Waveguiding properties of optical vortex solitons 
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We study the properties of linear and nonlinear waveguides induced by optical vortex solitons in 
both self-defocusing and self-focusing nonlinear media, for the case of a saturable nonlinear response. 
We demonstrate that the vortex-induced waveguides can guide both fundamental and first-order 
guided modes which together with the vortex may form, for larger amplitudes of the guide modes, 
different types of composite vortex-mode vector solitons. In the case of the self-focusing saturable 
medium, we demonstrate that a large-amplitude guided mode can stabilize the ring-like vortex 
structure which usually decays due to azimuthal modulational instability. 

PACS numbers: 



I. INTRODUCTION 

Vortices are fundamental localized objects which ap- 
pear in many branches of physics, ranging from the 
physics of liquid crystals to the dynamics of superflu- 
ids and Bose- Einstein condensates pj. Different types 
of vortices can be found and identified in optics; one of 
the simplest objects of this kind is a phase singularity in 
an optical wave front which is associated with a phase 
dislocation carried by a diffracting optical beam Q • 

When such an optical vortex propagates in self- 
defocusing nonlinear medium, the vortex core with a 
phase dislocation becomes self-trapped, and the result- 
ing stationary singular beam is known as an optical vor- 
tex soliton 3]. Such optical vortex solitons exist on 
a nonvanishing background wave, and they represent a 
two-dimensional generalization of the so-called dark soli- 
tons 0. They can be generated in experiment as self- 
trapped singularities of broad beams by using a holo- 
graphic phase mask, as was demonstrated for different 
types of nonlinear defocusing media 0,0,0- As a matter 
of fact, this type of optical vortex solitons demonstrates 
many common features with the vortices observed in su- 
perfluids and Bose-Einstein condensates @. 

In self-focusing nonlinear media, optical vortices can 
exist as ring-like optical beams with zero intensity at the 
center carrying a phase singularity J9|- However, due to 
the self-focusing nature of nonlinearity such ring-like vor- 
tex beams become unstable, and they are known to de- 
cay into several fundamental optical solitons flying off the 
main ring |10|. This effect was observed experimentally 
in different nonlinear systems, including the saturable 
Kerr-like nonlinear media, biased photorefractive crys- 
tals, and quadratic nonlinear media in the self-focusing 
regime (see details and references in Ref. |U). 

Waveguides induced by optical vortices in both lin- 
ear and nonlinear regimes are of a special interest be- 
cause this type of waveguides is robust and can be made 
rcconfigurable. Earlier numerical and experimental re- 
sults |ll|, Il2l [l3j indicate that many of the vortex waveg- 
uiding properties should be similar to those of planar 
dark solitons Q. Moreover, the vortex- induced waveg- 



uides can guide large-amplitude beams beyond the ap- 
plicability limits of the linear guided-wave theory, and, 
together with the guided beam, they can form a special 
type of the vortex-mode vector soliton or its various gen- 
eralizations HI El El E3. 

In this paper, we carry out a systematic analysis of the 
waveguiding properties of the vortex solitons and vortex- 
mode vector solitons in saturable nonlinear media, for 
both self-defocusing and self-focusing nonlinear media. 
We follow the earlier analysis of Ref. [13J and examine 
two major regimes of the vortex waveguiding. In the lin- 
ear regime, the guided beam is weak, and the analysis 
of the vortex waveguiding properties is possible by ap- 
proximate analytical methods, reducing the problem to a 
well-known analysis of the linear guided-wave theory (as 
an example, see Ref. E3)- The most interesting nonlin- 
ear regime corresponds to large intensities of the guided 
beam, and it gives rise to composite (or vector) solitons 
with a vortex component, that should be identified and 
analyzed numerically. 

We also analyze the vortex waveguiding and vortex- 
mode vector solitons in a self-focusing nonlinear medium 
where, as is well known, the vortex beam becomes self- 
trapped and it also displays strong azimuthal modula- 
tional instability. However, we demonstrate that a mu- 
tual incoherent coupling between the vortex waveguide 
and a large-amplitude guided mode it guides can provide 
a strong stabilizing mechanism for stable two-component 
vortex solitons to exist in such media. 

The paper is organized as follows. The next section 
(Sec. II) presents our model which takes into account 
both the incoherent coupling between the waveguide and 
the localized mode it guides, and the saturable nature of 
nonlinearity usually realized in experiment. The model 
is valid for both self-defocusing and self-focusing media. 
Section III is devoted to the study of the properties of 
optical vortex solitons and their guiding properties in 
self-defocusing media, whereas Sec. IV summarizes our 
results for the vortex waveguides and vector solitons in 
self- focusing nonlinear media. In particular, for the first 
time to our knowledge, we demonstrate that the vortex- 
mode vector solitons formed by the vortex waveguide to- 
gether with a large-amplitude mode it guides can show 
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the bistability property, when two types of different com- 
posite soliton structures exist for the same value of the 
soliton power. We also discuss, by means of direct nu- 
merical simulations, the stability properties of the vortex- 
mode vector solitons in a self-focusing nonlinear medium. 
In particular, we reveal an effective method of suppres- 
sion of the azimuthal instability through the stabilization 
effect that the guided mode exerts on the hosting vortex, 
the latter is known to be unstable by itself in focusing 
media. Finally, Sec. V concludes our paper. 



II. MODEL 

In order to study the guided modes of a vortex-induced 
waveguide in both self-focusing and self-defocusing me- 
dia, we consider the interaction of two mutually incoher- 
ent optical beams propagating in a nonlinear saturable 
medium. The equations of motion for two beams can be 
presented in the dimensionless form as follows, 



.du . i\u\ 2 + u\v\ 2 )u 

l— + Aj_U + T) /' ' , |2 I |2' 

oz 1 + a{\uY + \v\ z , 



dv . (lul 2 
i— + Aj_v + 7] 
oz 



fi\u\ 



r 



1 + a(\u\ 2 + \v\ 2 ) 



■0, 
0. 



(1) 



where u and v are the dimensionless amplitudes of two 
fields, the parameter a characterizes the nonlinearity sat- 
uration effect, and the incoherent mode interaction is de- 
scribed by the coupling parameter fi. The spatial coor- 
dinate z is the propagation direction of the beams, and 
stands for the transversal part of the Laplacian op- 
erator in the cylindrical coordinates r = (x 2 + y 2 ) 1 ^ 2 and 
<p = t&n (y/x). The sign parameter r\ = ±1 defines 
the type of the nonlinear medium under consideration, 
which is self-focusing, for r\ — +1, or self-defocusing, for 
77 = — 1, respectively 

The model Q provides a straightforward generaliza- 
tion to a number of important cases studied earlier. In 
particular, the limit a — > corresponds to the Kerr 
medium discussed, for example, in Refs. 0, 0, Il8j . 
whereas the saturable case at fj, = 1 corresponds to the 
incoherent beam interaction in photorefractive nonlinear 
media where different types of composite vector solitons 
have been predicted theoretically and observed in exper- 
iment 0,0,0,01- In particular, some examples of 
the two-component vector solitons composed of a vortex 
component and the localized modes it guides were pre- 
sented earlier by Haelterman and Sheppard [18( for the 
defocusing Kerr nonlinearity. 

We look for stationary solutions of the system ^ in the 
form of a radially symmetric single-charged vortex beam 
in the field u of the form 



u(r, 4>; z) — u(r)e 



(2) 



where 7 is the propagation constant of the vortex mode, 
which if considered normalized takes the values +1 or - 
1 in the focusing and defocusing cases respectively, and 
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FIG. 1: Examples of localized solutions for the vortex waveg- 
uide (u) guiding a fundamental- mode beam (v). The labels 
correspond to the points marked in Fig. [5] Two examples 
labelled by the same letter A constitute two different types 
of localized solutions found in this domain (mode bistability). 
The parameters are: A (fj, — 1.5, P = —1.1), B (p = 1.1, 
P = -0.85), C (p = 1.1, (3 = -1.0), D (fi = 2.0, f3 = -1.3), 
and E (p = 2.0, f3 = -1.05). 



consequently we have 7 = r\. For r — > 00, the amplitude 
u(r) approaches a constant value, for the self-defocusing 
medium, or zero, for the self-focusing case, respectively. 
We assume that the vortex guides the second beam given 
by the expression 



v(r, <f>; z) — v(r)e 



ilcji i{3z 



(3) 



where f3 is the propagation constant of the guided beam. 
In the expressions @, © the functions u(r) and v(r) are 
the radial envelopes of the interacting fields, and I is the 
angular momentum of the guided mode. 

In this way, the following set of stationary z- 
independent equations are obtained: 



. 1 (u 2 + av 2 )u 

-r]u + A r u -u + r\ - i _ t _ 9 i ^ = 0, 



I 2 

r 2 ' 



1 + cr(u 2 + V 2 ) 



(v 2 + [iu 2 )v 
-0v + A r v--v + v / \ ' = 0, 
r z 1 + <J{u + u ) 



(4) 



where A r = (l/r)d/dr(r d/dr). The model 0] describes 
different types of localized solutions for the vortex and 
the localized modes it guides, and it has been further 
analyzed by the numerical relaxation methods to classify 
its localized solutions for both defocusing and focusing 
cases, respectively. 
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In the defocusing case when rj = — 1, we have u — > uq 
and v — > for r — > oo. Using this asymptotic conditions, 
we obtain 



1 



1 + CTUq 



u = 0, 



(5) 



and find the result for the vortex background amplitude 
wo in the form, uo = (1 — cr) -1 ' 2 . 

Using the second equation for u, and neglecting the 
terms which most quickly tend to zero when r — > oo, we 
obtain: 



-f3v- 



cPv 
dr 2 



ml 



1 



0. 



(6) 



In order to have the exponentially decaying function v 
(for r — > oo), the following condition should be satisfied, 



-(3- 



1 + CTWq 



< 0. 



that gives as the result 

(3 + ^>0 



(3 > -fi, 



(7) 



(8) 



which defines the existence domain threshold as the line 
—f3 = fi. This result is valid for both fundamental and 
first-order guided modes. 



III. VORTEX GUIDED MODES 

First, we consider the case of self-defocusing nonlin- 
earity when rj = — 1. In this case, the vortex amplitude 
u{r) approaches a finite-amplitude asymptotic value uq 
for r — > oo, where uo can be found from the asymptotic 
analysis, Uo = (1 — cr) -1 / 2 . With the help of numeri- 
cal relaxation methods, we find various types of localized 
modes that describe a vortex and the localized guided 
mode with the amplitude v(r) satisfying the condition 
v(r) — > for r — > oo. This condition requires that (3 < 0, 
since the field v(r) is localized due to the presence of the 
vortex waveguide only, and it describes a guided mode. 

Figure ^ shows several characteristic solutions of this 
type for particular values of fj, and a = 0.5, whereas the 
existence domain for the solutions with the fundamental 
guided mode (i.e. when I = 0) is shown in Fig. [21 For 
the values of the propagation constant (3 close to the 
cutoff (f3 c ~ —0.79 at fi = 1), the second field presents 
a low-amplitude guided mode that can be analyzed by 
means of the linear guided wave theory. However, for 
values of (3 close to the threshold, the amplitude of the 
guided mode becomes comparable to that of the vortex 
and its width grows, as shown for the case E in Fig. ^ 
In this latter case, the vortex is deformed as well, so the 
resulting structure describes a kind of a two-component 
vector soliton with the vortex-like mode. 

In Fig. we show the existence domain for the first 
vortex-like (1 = 1) mode of the induced vortex waveguide. 
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FIG. 2: Existence domain of the fundamental mode in the 
vortex-induced waveguides on the plane (fj,, (3), for a = 0.5. 
Marked points correspond to the solutions presented above in 
Fig. A triangle region with the point A corresponds to the 
existence domain for bistable composite solitons. 
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FIG. 3: Top: existence domain for the first-order mode (I — 1) 
guided by the vortex-induced waveguide, show on the plane 
(fj., P). Below: examples of the two-component localized so- 
lutions for a = 0.5, fj, = 2, and (3 = -1.05 (A), (3 = -1.2 (B), 
and P = — 1.9 (C), respectively. 



Similar to the case of the fundamental guided mode, the 
linear theory is applicable near the mode cut-off, as is 
shown for the mode C, where the mode grows and it 
deformes the vortex waveguide substantially, see the case 
A near the upper boundary of the existence domain in 
Fig. El (top). 



IV. SELF-TRAPPED VORTEX BEAMS 

A. Vortex-mode vector solitons 

In the case of the self-focusing nonlinear medium, we 
put r) = +1, and consider the self-trapped vortex beams 
with the asymptotic behavior u(r) — > as r — > oo. Be- 
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FIG. 4: Examples of solutions for the fundamental mode of a 
waveguide induced by a bright-like vortex. Labels correspond 
to the points marked in FIG. 4, using a subscript to name 
the two solutions found for the points B and C. Values: A 
(At = 1-10, (3 = 1.38), B (ai = 1.10, P = 1.44), C(p= 1.17, f3 = 
1.38), D (At = 1.17,/?= 1.47). 



sides that, it is required that v(r) — > and (3 > 0, so that 
v(r) describes a localized mode. At the origin (r = 0), 
the boundary condition for the vortex is u = and, if 
we seek the solutions corresponding to the fundamental 
guided mode (1 = 0), we have v = vq (or dv/dr = 0) for 
the field v(r) at r = 0. 

In Fig. we show some examples of the two- 
component localized solutions which describes a funda- 
mental (no nodes) beam guided by the self-trapped vor- 
tex in the form of a vortex-mode soliton. The existence 
domain for this kind of solutions is shown in Fig. \5\ on 
the parameter plane (a*, /?), calculated numerically for 
a = 0.5. Close to the cutoff, the guided mode has a 
small amplitude (see the cases A and D2), and the vor- 
tex is only weakly distorted. However, when the prop- 
agation constant (3 is close to the upper threshold, the 
field v(r) grows in the amplitude affecting strongly the 
vortex mode. 

In order to describe the bistable vector solitons, in 
Fig. we plot the bifurcation diagram of the two- 
component localized solutions, for the partial and to- 
tal beam powers. The power of the composite vortex- 
mode solitons Ptotai originates at the bifurcation point 
0\ where the mode v is small and it can be described 
by the linear theory. For larger value of /3, the curve 
bends, and then it merges with the other partial power 
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FIG. 5: Existence domain in the plane (a*,/?) of the fundamen- 
tal mode guided by the self-trapped vortex for a = 0.5. Inset 
shows a zoomed region where two different types of solutions 
exist (bistability zone). 
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FIG. 6: Bifurcation diagram for the one-and two-component 
solitons composed by the bright vortex inducing a waveguide 
together with its fundamental mode for a* = 1.1. Thin lines 
represent the powers P u and P v of each of the scalar solitons; 
the thick line is the total power Ptotai of the composite soliton. 
The points Oi and O2 are the bifurcation points where the 
vector soliton emerges, and the points Bi and B2 correspond 
to the examples of the bistable solutions for ft — 1.44 shown 
in Fig. H 



curve P v at the bifurcation point O2 (see Fig- EJ). The 
bistable solutions Bi and B2 correspond to a single value 
of the propagation constant (3 in the bistability domain. 
Importantly, two solutions have different stability prop- 
erties, and only one of the solutions is stable, as shown 
in Fig. □ 



B. Vortex stabilization 

The incoherent interaction between the vortex and the 
localized mode it guides has the character of attraction, 
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FIG. 7: Examples of the vortex evolution in the bistability 
domain. Shown are the field intensity profiles as gray-scale 
images at several propagation distances. Top row: the com- 
ponents of the vector soliton labelled as Bi in Fig.HJ Bottom 
row: the components of the vector soliton labelled as B2 in 
Fig.© 



and it may provide an effective physical mechanism for 
stabilizing the vortex beam in a self-focusing nonlinear 
medium. Indeed, it is well-known that the vortex beam 
becomes unstable in a self-focusing nonlinear medium 
due to the effect of the azimuthal modulational insta- 
bility. In this case, the vortex splits into the fundamen- 
tal beams that fly off the main vortex ring [lOj . On the 
other hand, the bright solitons are known to be stable in 
such media. We expect that a mutual attraction of the 
components in a two-component system would lead to 
a counter-balance of the vortex instability by the bright 
component if the amplitude of the latter is large enough. 

To carry out this study, we consider a two-component 
composite structure consisting of a vortex beam together 
with the fundamental mode of the waveguide it guides, 
described by Eqs. ^ ( a t V = 1 ancl M = 1-0)- To study 
the mode stability, we propagate the stationary soliton 
solutions numerically. 

In Fig. |S1 we compare the vortex stability for the scalar 
and two-component systems. In the top row, we show 
the propagation of a vortex alone in the scalar model; 
the vortex breaks into two solitons which fly away after 
some propagation distance. In the bottom row, we show 
the propagation of two coupled components (the vortex 
and bright mode it guides). Due to a strong coupling be- 
tween the modes, the propagation of the vortex is stabi- 
lized and its decay can be delayed dramatically, as shown 
in Fig. flower row), or become completely stable, simi- 
lar to the case shown in Fig.|7Jlow row). We confirm this 
stabilization mechanism by performing a similar study 
for a Gaussian input beam of the bright component in- 
stead of the exact stationary state, as is easier to realize 
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FIG. 8: Vortex stabilization in a two-component system. 
Shown are the field intensity profiles as gray-scale images at 
several propagation distances. Top row: propagation of a vor- 
tex soliton by itself (v — 0). Bottom row: both the vortex 
beam and fundamental mode (approximated by a Gaussian 
profile) propagate together. Parameters are: a = 0.5, /x = 1, 
and p = 1.45. 

in experiment. 



V. CONCLUSIONS 

We have studied the properties of the vortex- 
induced optical waveguides in both self-focusing and self- 
defocusing saturable nonlinear media. We have calcu- 
lated numerically the existence domains of the compos- 
ite solitons created by an optical vortex and the local- 
ized mode it guides, and have analyzed the properties 
of the vortex waveguides in both linear and nonlinear 
regimes. In the case of self-defocusing nonlinear media, 
we have identified the regimes when the vortex soliton 
can guide both the fundamental and first-order guided 
modes creating an effective multi-mode waveguide. In 
the case of self-focusing nonlinear media, we have stud- 
ied the effect of the strong coupling between the vortex 
waveguide and the guided mode on the vortex stability. 
In particular, we have described a novel mechanism of 
stabilizing the vortex azimuthal modulational instability 
by a large-amplitude guide mode, and we have revealed 
the existence of optical bistability for the vortex-mode 
composite solitons. We believe similar results can be ob- 
tained for vortices in other types of nonlinear models, and 
they can be useful for other fields such as the physics of 
the Bose-Einstein condensates of ultra-cold atoms. 
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